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ABSTRACT 



The optimal time- sequential distribution of supporting 
fire against enemy ground units in combat against attacking 
friendly units is studied. Lanchester- type models of war- 
fare are combined with optimal control theory in this in- 
vestigation. The optimal time-sequential fire-support 
policy is characterized for a specific problem. Although 
complete details for the determination of the optimal policy 
are not given, it is conjectured, on the basis of the 
theorems which were proved, that for this problem the optimal 
policy is to always concentrate all supporting fire on the 
same enemy unit (until supporting fire must be lifted). 
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I . INTRODUCTION 



The problem o£ fire support allocation is a problem of 
interest to the military tactician and planner. In this 
paper, we will study the problem of artillery fire support 
allocation against several enemy ground units. The problem 
is to determine an optimal time-sequential policy for dis- 
tributing available fire support without wasting it. By 
wasting of fire support we mean overkilling, i.e., the 
destruction of enemy forces that does not contribute to the 
attainment of the friendly objective. Overkill would cor- 
respond to a state variable which represents the strength 
of a ground unit becoming negative. We do not want to over- 
kill because it will be simply wasting the effort. 

To determine a "good" allocation of fire support, one 
must consider the dynamic nature of combat in which opposing 
forces trade casualties over a period of time. Lanchester 
type models of warfare (see [15]) have been developed to 
provide insights into the dynamics of combat. Although com- 
bat attrition certainly is a random process, such a formula- 
tion is not mathematically tractable as far as analytic v^^ork 
is concerned. Therefore, we will consider deterministic 
Lanchester- type models in our work here. These are differen- 
tial equation models of the combat attrition process. With 
such a dynamical system model and considering that fire sup- 
port decisions are made over a period of time, one is led 
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to an optimal-control/di££erential-game £ormulation £or the 
determination o£ optimal time-sequential £ire support 
allocation. 

Such time-sequential combat optimization problems may 
be cast as either optimal control problems or di££erential 
game, depending on whether or not both sides are modelled 
as rational decision makers. In the exploratory work at 
hand, we will consider the optimization o£ combat decision 
by only one o£ the combatants. This leads to consideration 
o£ an optimal control problem. Such a problem has been 
re£erred to as a Lanches ter- type optimal control problem by 
Taylor (see [14] ) . 

In a deterministic combat optimization problem, we usually 
encounter either a static or dynamic situation. I£ the sit- 
uation. I£ the situation is static, we may determine the 
optimal policy using the classical techniques £rom advanced 
calculus £or problems that are either unconstrained or sub- 
ject to equality constraints. I£ there are inequality con- 
straints in the optimization problem, the methods o£ nonlinear 
programming may be used. I£ the combat situation is dynamic, 
we may apply the methods o£ the classic calculus o£ variations 
£or problems that are either unconstrained or subject to 
equality constraints (see [4]). However, i£ inequality con- 
straints arc involved, it may be necessary to apply the 
recent results o£ optimal control theory. 

The development of solutions to Lanchcstcr- type optimal 
control problems is not a routine matter due to several tech- 
nical difficulties. Due to physical considerations force 
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levels cannot become negative, and we must introduce state- 
variable inequality constraints (SVIC's), i.e., no state 
variable (force level) can be allowed to be negative. All 
Lanchester- type dynamic tactical allocation problems have 
SVIC's. The maximum principle (in its original form) is 
inadequate to solve a problem with SVIC's, since an ex- 
tremal may contain a subarc which lies on the boundary of 
the state space. By an extremal we mean any trajectory on 
which the necessary conditions of optimality are satisfied 
for all t. IVhen SVIC's are present in a problem, the 
necessary conditions (i.e., Pontryagin's maximum principle) 
need to be modified (see [7]). To simplify considerations 
for the problem studied in this thesis, we will assume that 
conditions are such that all state variables are always 
non-negative. Thus SVIC's are not considered further in 
this thesis. 

Another difficulty for this type of problem is the pos- 
sible presence of a singular subarc. Consider a control 
problem in which, for example, a single control variable u 
appears linearly in the Hamiltonian H. The maximum prin- 
ciple fails to determine the optimal control when 3H/8u = 0 
for a finite interval of time. The corresponding trajectory 
is called a singular subarc . Use of the terminology singu- 
lar is motivated by the fact that in such problems the matrix 
of second partial derivatives of the Hamiltonian with respect 
to the control variables is a singular matrix. We will call 
an extremal singular if it contains a singular subarc. 
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We will study the problem at hand in two phases. In the first phase 
we will consider nonsingular extremals, while in the second phase we 
will consider singular extremals. One result of the study for this 
problem is that there is no nice algorithm to determine the optimal 
control for our problem. In general, combat optimization problems 
dealing with multiple units are extremely complicated, especially if 
there are some singular subarcs. 
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II. ARTILLERY FIRE SUPPORT ALLOCATION MODEL 



Consider combat between two ground forces Red and Blue. 
Each side is composed of two infantry units. Blue infantry 
units are being supported by artillery fire. 

We want to find the best policy for allocating Blue's 
artillery fire support against Red to achieve the maximum 
effectiveness. The measure of effectiveness is taken to be 
the ratio of Blue infantry forces to Red at the given time. 
The problem is to determine an optimal policy by which the 
ratio of two infantry forces at the end of the battle can 
be maximized. 




•A' 
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The Blue forces will be denoted as and X^2 (with 



corresponding force levels and x^2Ct)) and the Red 

forces as X22 and X22 (with corresponding force levels 
X2j^(t) and X22 Ct) ) . 

The following assumptions will be made for developing 
the model. 

1 . Each infantry unit is only in contact with his op- 
posing unit. 

2 . The Blue infantry units are moving to contact with 
the Red units and due to their movement cause negligible 



attrition to the opposing Red forces. Furthermore, the 
Red defensive positions are prepared and relatively invul- 
nerable to small arms fire during this approach phase. 

5 . The defensive Red units cause attrition to the op- 
posing Blue units according to a "square law" attrition 
process . 

4 . The Blue artillery delivers area fire against the 
Red infantry. 

With the above assumptions we may develop the following 
combat optimization model: 



maximize 



X (T) X-.2CT) 



subject to: = '^ 1 ^ 21 ^^^ 



^12 Ct) 




X22(t) = -C2(l - u(t))x22(t) 



(la) 
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I 



where 



x..rtl: force level at time t 

ij ' 

X- • ft) : the time derivatives of state variables x..(t) 
T: time at which the battle terminates 
uft) : control variable 
a« , c*: attrition rate coefficients 

: utilities assigned to ratio of survivors. 

Since the criterion functional contains force ratios at 
the terminal time, we may formulate a problem with a state 
space of lower dimension by defining new state variables as 

ill = X = X . 

^21 ^ ^22 ^ 

Then we can show that (la) is equivalent to the following 
optimal control problem; 

maximize J = kj^x^(T) + k 2 X 2 (T) 

subject to x^ = -a^ + c^ux^ 

X 2 = -Si2 (lb) 

where all constant coefficients remain same as fla) . 

Let us assume that we are given the initial conditions 

x^(O) = x° , X2(0) = x°. 

Let us assume also that we are given a specified time which 
is called the duration of the combat. We will assume that 
the battle will be terminated at time T under one of the 
following conditions: 
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Condition 1; x^(T) = X 2 (T) = 0, for some T<T 



1 * 



Condition 2: T=T 



III. DEVELOPMENT OF NECESSARY CONDITIONS OF OPTIMALITY 



To develop the necessary conditions for optimality, let 
us define the Hamiltonian for the problem (lb) , 

H(t,x,p,u) = p^(-a^ + c^uxj^) + P 2 (-a 2 +C 2 (l-u)x 2 ) ; (2) 



Pj^ and P 2 are the adjoint variables, corresponding to the 
state variables x^ and X 2 , respectively and, satisfying the 
differential equations: 



‘^Pl _ 3H _ 

“at " ' ^ “ -c^up^ 



dp 2 
dt 




= -C2(l-u)p2. 



(3) 



By Pontryagin's maximum principle, we can show that the 
nonsingular optimal control will be obtained by maximizing 
the Hamiltonian over the control variable. Hence we may 
observe that the extremal control is then given by 

,0 if S < 0 

u(t) = ] 1 if S > 0 

V indeterminate if S = 0, 

where S = ‘ ^2P2^2‘ 

We will call u the extremal control function and S the switch- 
ing function. The extremal control u so defined may furnish 
the global maximum, a local maximum, or only a stationary 
values for the criterion function. 
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Since the Hamiltonian H is a linear function of the con- 
trol variable u, if the switching function S becomes iden- 
tically zero for some finite interval of time, the maximum 
principle fails to provide any information about the optimal 
control u. In this case all the admisible values of u 

maximize the Hamiltonian; therefore theory of singular ex- 

3 H 

tremals is required to solve this problem. If S = is to 
be identically zero for some finite interval of time, then 
all its derivatives with respect to the time must also van- 
ish. That is, the trajectory remains on a singular subarc 
if and only if 



m ^ r lii 'i = r M 1 
3u 3T ^ 3u dt2 '■ 3u ^ 



• • • 



(5) 



From these conditions, we will have 



0 S c^p^x^ - c.^p.^x 



l^^l^l ^-2^2 2 



® ^ 3u ^ ^ ‘^l^l^l ^2*^2P2 



0 - ^ ^ - S - a^c^p-,u - a^c^p^(l-u) 



‘l^l^l' 



From (6) - (8) , we have 



( 6 ) 

(7) 

( 8 ) 



u = 




s , 



(9) 



defining s, the singular control. 

The generalized Legendre -Clebsch condition is necessary 
for singular subarc to yield a maximum return (see [9]). 
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The condition will be given 
the problem at hand. 






( ^ )] 
^ 3u ^ 



< 0 . 



For 



9 r d^ 

9u *■ dt^ 



Mi 

9u 



) ] = 



a^c^Pl + a2C^p2 > 



0 , 



since p^ > 0. Therefore, the Legendre-Clebsch condition is 
satisfied. 

For the singular subarc S = S = 0; it follows that 



x^Ct) 



( 10 ) 



We will refer to the surfaces S defined by (10) > and S , 

X p 

defined by (7) as singular surfaces. Singular arcs lie in 
these. The surface defined by 



S = - C 2 P 2 X 2 = 0 



is called the switching surface. 

The surfaces or curves S and S,, can be drawn in the p 

p X ^ 

and X space, respectively, but the sxvitching surface S in- 
volves both p and x and we cannot visualize it easily. 
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IV. SYNTHESIS OF EXTREMAL CONTROL 



A. SYNTHESIZING FUNCTION 

Let X be the Euclidean space on which the state variables 
X and their associated adjoint variables p are defined, and 
U be the Euclidean space for control variables. Suppose we 
can find a real valued function g: X-»-U such that x = f(x,g(x,p)j 

determines all the extremal trajectories from initial to ter- 
minal time. Then the optimal control function may be explicit- 
ly determined by the function of state variables and adjoint 
variables such that the optimal control at each time depends 
only on the points of the space at which the state variables 
and the adjoint variables are located at the given time. 

The function g(x,p) will be called a synthesizing function 
and the explicit determination of such a function is called 
the synthesis of the extremal control, and the explicit de- 
termination of the time history of the extremal trajectory 
is called synthesis of extremal trajectory. 

Unfortunately, an existence theorem for the synthesizing 
function has not been proved. Furthermore, there is no general 
algorithm for the problem of finding such a function. If the 
set of terminal values for the state variables is not explicit- 
ly specified, it may be extremely difficult to find a syn- 
thesizing function. The determination of the synthesizing 
function must be done on a case-by-case basis, since there is 
no general algorithm. In the problem at liand the terminal 
state is not specified. However, we know that the control 
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depends upon the state variable x and the adjoint variable 
p, and it will be determined by the maximum principle when 
(x,p) does not belong to the singular surface. 

The difficulty is that x is only known initially and p 
is known only at the problem's end until we have managed to 
match X and p with the maximum principle to get a solution. 

However we may analyze the behavior of the switching 
function, which is a type of synthesizing function. Although 
it involves both state and adjoint variables we can infer 
its behavior as a function of time, along an extremal, to 
establish several results. Further the trajectory lies on 
singular subarc for some finite interval of time, we may 
consider the singular surfaces and as a synthesizing 
function. 

B. SYNTHESIS OF CONTROL FOR A TIME INTERVAL DURING WHICH 

THE CONTROL IS CONSTANT 

For synthesizing extremal trajectory, let us consider 
the constant control at first. Suppose we choose a constant 
optimal control from the beginning to the end of the battle. 
Then from the differential equation (lb) , we will have the 
following solutions. 

For u = 1 ; 

x^(t) = max ^ )e ^ , 0 ] 

X2(t) = max [ x^ - 0 ] (12a) 
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PlCt) = k^e 



c^CT-t) 



PzCt) 

S(t) 

S(t) 

S(t) 



= k. 



Ci(T-t) 



-a^c^k^e 






c^(T-t) 
c^(T-t) 



* 



For u = 0 ; 



= max {x^ - a^t, 0} 



X- (t) = max { (x 



3./^ O ^ 3^ 

0 ^ \ ^ ^ _L ^ 



)e " + , 0} 



PlCt) 



= k. 



P2(t) = k2e 



C2(T-t) 



set) = - CjkjXje 



S(t) = + a2C2k2e 

S(t) = -a2C2k2e 



C2(T-t) 

CjCT-t) 



For u = 3 = 



Cl + C2 



3 3 

Xi(t) = max {(x° 1 “I » 



X2(t) = max { (x° | )e“^ + — 



a 



, 0 } 



Pl(t) = kie“^^‘^^ 
P2(t) = k2c“(T-t) 



(12b) 

(12c) 

(12d) 

(12e) 



(13a) 

(13b) 

(13c) 

(13d) 

(13e) 



(14a) 



(14b) 
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(14c) 



S(t) = ic^k^x^ - C2k2X2)e“^’^'^^ 



S(t) = C-a^c^kj^ - a2C2k2>e“^’^‘^^ 



(14d) 




(ajcJkj - a^c^k^) 



C14e) 



where 



a 




We will use the above results to determine the optimal 
policy for more general cases. 

C. DEVELOPMENT OF THEOREMS ON THE CHARACTERIZATION OF AN 
OPTIMAL POLICY 

By analyzing the behavior of switching function, we may 
establish the following theorems. 

THEOREM 1 

Let E be a maximizing extremal which contains no singular 
subarc for the optimal control problem given by (lb) . Then 
there is at most one switching in E. 

Proof : 

By hypothesis, since E does not contain a singular sub- 
arc, the optimal control can be determined by Pontryagin’s 
maximum principle, except possibly at one point, as we shall 
see. We have got the switching function from the Hamiltonian 
by means of the maximum principle. 



S(t) > 
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Then its derivative is 



S(tJ = ajC2P2 - a^CjPj 



and 



S(t) = a^c^p^u - a2C^P2(l-u). 



It follows that 



*S*Ct) > 0 if u = 1 
S(t) < 0 if u = 0. 



Since the trajectory does not lie on singular subarc, we will 



have either u= loru=0att=0. Let us consider the case 
u=latt=0;a similar discussion will hold if u = 0 ini- 
tially. Since u = 1 at t = 0, S(0) > 0. 

1. Let us consider first the cases in which S(0j > 0. 
There are two sub-cases, accordingly as S(0j > 0 or S(0) < 0. 

Case 1 : S(0) > 0 and S(0) > 0. 

Since S(0) > 0, we must choose u = 1 at the beginning of 



fore, Pj^(t) is decreasing function of t and ?2 (t) is constant. 
Hence S(0) > 0 implies S(t) > 0 and S(t) > 0, for all t > 0. 
Therefore S(t) is a monotone increasing function, and conse- 
quently we will have S(t) > 0 for te[0,T]. Hence no switch- 
ing occurs. 

Case 2 : S(0j > 0 and S(0) < 0. 

From equation (12c) , S must have the form 



the game. From (12b), p^(t) = k^e 




, P 2 (i) = ^ 2 . There- 



- c. t 

D(t) = -Ae ^ + B 



( 15 ) 
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in which A and B are positive, so long as u = 1. But if S 
becomes negative u must become 0 and S will have a different 
form. 

Let T be the value of t for which D = 0. Consider also 
S: S may become zero at some value t^ < t or S(t) may be 
greater than or equal to zero. There are three subcases. 




It may be that S > 0 for all t < r. In this case S > 0 
for all t, since S > 0 when u = 1. Hence no switching occurs. 
Subcase ii 

It may be that S becomes zero at some time t^ < t and 
S(t^) < 0. Then there is a switching at t = t^. From the 
facts that S(tj^) = 0,S(t^) < 0 and S < 0 when u = 0 it fol- 
lows that t^ is a switching point, and that S < 0 for all 

t > tj. 

Subcase iii 

It may be that t^ in case ii coincides with t, so that 
S(tj^) = 0 and S(t^) = 0 for some t^^. < T. Then the maximum 
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principle does not determine the control at time t^. If we 

• • 

choose u = 1 at t. = t^, when S(t^) >0. It follows that 
S(t), S(t) are positive for t > t^^ and there is no switching. 
If we choose u = 0 at t = tj^, then by the above arguments 
S(t) < 0 for t > t^ and there is one switching. 

2. It may be that S(0) = 0. If S(0) > 0 then S(t) is 
always positive. If S(0) < 0 then S(t) < 0 for t > 0 and 
u = 0. Finally if S(0) = 0 then we may choose u(0) arbi- 
trarily. If we choose u = 1 initially then S > 0 for t > 0 
by argument given in case 2-iii and there is no switching. 
Similarly, if we choose u(0) = 0, S remains negative. 

This takes care of all extremals which do not have a 
singular arc; a nonsingular extremal has at most one switch- 
ing point. 

COROLLARY ; 

If E contains a singular subarc, then there are at most 
two switchings. 

Proof ; 

Suppose S(0) > 0 and there is tj^ > 0 such that S(tj^) = 
S(t^) = 0. We will show that there are at most two switch- 
ings for this case. The case for S(0) < 0 will be symmetrical 
to the case of S(0) > 0. Since we will choose the singular 
control at t^, we may assume that the trajectory will remain 
on the singular subarc until some time t 2 » If t 2 > T, we 
will have only one switching at t^ from u = 1 to u = 

Suppose t^ < T. Assume that at some time t^ we choose u = 1 
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when we are on the singular subarc. Then, by the argument 
given in 2-iii, S will become positive and remain positive 
for all t > t 2 « Similarly if we choose u = 0, S will be- 
come negative and remain negative. Thus we will have at 
most two switchings, one at t^, say, when the extremal enters 
the singular arc and one at t 2 when it leaves it. 

THEOREM 2 

Consider a path E^q defined by the constant control 
u = 1, 0 < t < t^, and u = 0 , t^ < t < T , so that E^^ has 
one switching and no singular subarc. Let be the as- 
sociated criterion function. If > a^^/c^ and X 2 > & 2 T, 
then there is no maximizing extremal of the form E^q. By 
symmetry, there cannot be one of the form E^^ in which 
^2 - ^2^*"2 ^1 - ^1*^* That is, a properly chosen constant 

control is always better than a mixed policy of this type. 



Proof : 

Suppose we choose the extremal control u = 1 for te[0,tj^], 
and u = 0 for te[t^,T]. Then from (12a) and (13a), we have 

x^(t^) = max UxJ - )e ^ 

^2^^1^ = max { X 2 - a 2 t^) , 0). 



For te [t^,T] , 



x^(t) = max (x^(t) - 

^1 

= max {x? - — 

1 “"i 



a (t- 1^) , 0 } 



*" 1^1 

)e ^ ^ ^ 

^1 



(t-t^) ,0} 
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i 



^2 ^2 

x^Ct) = max{(x°(t^) - ^ )e ^ + — , 0 

2 2 

3-y a« 

= max{ (X2 - ^2^1 " 3 ^)® +3^>0}' 

2 2 

Suppose x^(t), X 2 (t) > 0 for all te[0,T]; then the payoff 
will be 
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3--I C-| t > a« 

ki[x° - ^ )e 1 ^ ^ - a^CT - x^)] + 



(16) 



kz [ (x° 



° ■ ^2^1 



^2 ^ ^2 



)e 



^2 



a- c^ti a, c,CT-t,) 

‘^10 “ ®iki[xj ■ ^ )® * ^ ' ^2^2^ (^Z'^2^1^ 

E-| c -1 1 -1 C/,(T"t-|) 

^10 = =i''lC='l - F7 * "2''2® ‘ 1 

a± x° 

Since x? > — and T < — , we always have JV^ > 0. Therefore, 

-L ” “ 3^2 

we will have a maximum payoff at the boundary value of t^, 
i.e. , t^ = 0 or t^ = T. 

If t^ = 0, 

nax . J„ - kjCxJ-a^T) * k2[x^ - ^ ) e ^ ] 

If = T, 



max Jjj = Jj = kj[x” - ^ )e ‘ ^ ' *2'^) • 



In other words, we cannot improve the payoff by switching the 
control. We will have the same result for Jg^j for which we 
switch the control from u = 0 to u - 1. We may summarize 
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this in the following corollary: 



COROLLARY : 

All nonsingular controls which maximize the criterion 
function for this problem have the form u e l or else u e 0. 
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V. DETERMINATION OF THE OPTIMAL CONTROL 



A. SUFFICIENT CONDITIONS FOR OPTIMALITY 

As is well known the maximum principle only provides 
necessary conditions of optimality. Second order sufficient 
conditions (see [1]) are impractical to apply to the problem 
at hand. Mangasarian ' s sufficient conditions (see [10]), 
unfortunately do not apply to the problem at hand. 

The following is an adoptation of Mangasarian ' s Theorem 
to a form more closely related to the specific problem under 
study. Consider the following optimal control problem. 

•'T 

FCt,x,u)dt + G(x^,t^) 
o 



Maximize 

U 



subject to X = fCt,x,u), 

gCt,x,u) > 0, 

the initial conditions C. fx ) >0, • 

the terminal conditions C 2 (x£) > 0, 

where x^ = x(0), X£ = x(T) for specified T. 

THEOREM : 

Assume that F[t,x,u) and each component of fCt,x,u) and 
g(t,x,u) are differentiable and concave in the variable (x,u) 
for te[0,T], and also each component of and C 2 (x£) is 

differentiable and concave in x^ and X£ respectively. Then 
a control u(t) yields a maximum for problem (17) if there 
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exist u(t) , x(t) , p(t), q(t), r(t) , s(t) such that; 

H(t,x,u,p,q,r,s) = F(t,x,u) + p^£(t,x,u) + q^g(t,x,u) 



+ rC^(x^) + sC2(x^) 




= F^(t,x,u) + p^£^ 
= F^(t,x,u) + p"^£^ 
(Xq,x£) 4- 7c^Cx^) 

(XqjX^) + sC2(x£) 

r > 0, 
rC- (x ) = 
s > 0 
sC2(x^) = 
q > 0 
qg(t,x,u) 



(t,x,u) + 
(t,x,u) + 

+ p(0) = 0 

- P(T) = 0 



0 



0 



= 0 



-T 

q g 



-T 

q g 



X 



u 



(t,x,u) 
(t ,x,u) 



+ pCt) = 0 
= 0 



p > 0 . 



However, £or the problem at hand (lb) it is readily seen 
that the right-hand sides o£ the di££erential equation are 
not concave £unctions o£ the state and control variables so 
that the above theorem cannot be applied. We will return to 
discuss our problem. 
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B. DETERMINATION OF OPTIMAL CONTROL IN TERMS OF INITIAL 
CONDITIONS AND TERMINAL CONDITIONS 

In general we cannot determine the optimal control with- 
out some computation and comparison. However, in some cases 
we can establish its form easily. In this section two theo- 
rems are proved. The first one establishes a sufficient 
condition for some nonsingular arcs to be optimal. The 
second theorem establishes a necessary condition: An opti- 

mizing arc never has a terminal singular subarc. 

THEOREM 3 

. ^1 

If E does not contain a singular subarc, and if x. > — 

clcx^ 1 

and 0 < T < min{ In constant control 

^2 ^1 

u = 1 is the optimal control. By symmetry, if X 2 > — and 

X ® c Ic X ® 2 

0 < T < min{^ ^ iln constant control 

3.^ f 

u = 0 is the optimal control. 

Proof : 

^1 

Let E be a nonsingular extremal for which > 7 r~ and 
x° • 1 

— > T. Then, by theorem 2, we know that an optimal control 

^2 

is constant; either u e 0 or u = 1. Now assume S(t) < 0 
for all te [0 ,T] . 

c.T 

From Cl3c) , S(0) = -C 2 k°e + c^k^x® < 0. Hence 

c k X ® 

X > 1 r . Therefore, by contrapositive argument, 

if T < (l/c 2 )S-n(c^k^xJ/c 2 k 2 Xp , then S(0) > 0. Thus u = 0 
is not an optimal control, consequently the optimal control 
is u = 1. A similar argument establishes the symmetrical 
result . 
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THEOREM 4 



A maximizing control does not have a singular subarc 
ending at t = T. 



Proof : 

Assume that we have a path which has singular control at 
the end, say, at least in an interval t^^ < t < T. Let us 
consider At = T - t^ to be small. We want to compare this 
path with paths which have other controls in the interval 
(tj,T) . In particular we want to compare the three trajectories 
which have the controls u= s, u= 1, oru= 0, in the final 
interval (t^,T). 

Consider the Taylor expansions of J at t = t^. 



dJ 



d^J 






at " I ( 



). ‘At^ + O(At^) 



dJ, T d^J^ 



dJ, 



Jq = * ( 



dt t 



1 

-At + I C \ *At2 + O(At') (18) 



Since we assume that the trajectory lies on the singular sur- 
face at the end of the battle, to have singular control on an 
admissible extremal at time T, it is necessary that 



S(T) = S(T) = 0, 



From (6) and (7) , we have 



a^Ciki = , i.e., k^ - j^kj 
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Also 




^1 ^2 

in the interval of singular control. Then at time t^, this 
condition must be satisfied. 

From (10) , (12a) , (13a) and (14a) , we may write the 
payoff 



a. a(T-t^) 

+ rxi i le ^ 



) [ — + (x^ )e ■" ] 

s l^ao '^'•a^l •' 



where = x^^(t^) > 0 , x^ = X 2 (t^) > 0. 



Similarly, 



a. c^ 

[ — + — xj 
1 1 L C2 1 



^1^1 1 ^1 ^iCT-t^ 

(T-t.^) + (ij - ^ )e 1 



Jo = >^1 [ 



a.c. . c. ^2^"^”^!^ 

+ XT - a^(T - t.) + ( — Xt - ' 2 - )e ^ '^ ] 

c| 11^ 1^ ^ C 2 1 c| J 



Their derivatives with respect to t^ are; 

a^ a(T-t.) 
j; = -kjCj (I* - ji )e 

a. a(T-t.) 

Js ' ^ )e 

^1^1 t ^1 

Ji = ->'1 [ - 4^ * ' 

a, c,(T-t,) 

Ji = V!(='i ■ 

■>0 ‘ ->'1 [ -“1 " CCiX} - )e " J ] 
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a- c-(T-t-i) 

j" = k. c c fi* - — 1e ^ ^ 

•^0 1 . Co 



Hence , 



Vz 



1% < = 



lim J\' = k-c?(x|^ - — ) 

^ .-r 1 1 1^ 1 



t^->-T 



lim J" = 
t^-^T ^ 






Since ~ ^'^O^T’ have 



= Irm 



3^1 



Jl - Jg = J ^’^1 • 



3 _ 



Cj * Cj 



> 0 



Jq - Js = lfj() - J's’CT - ti) - 0(T - tp 



k, c-Coxf 



3 , 1 ^ 1 ^ 2-1 , 0 ^ 

1 ■ ''2 



c, + c. 



Therefore, J- > J , and J > J , at least if T - t-, is suf- 
ficiently small. 

Thus we see that if a trajectory ends with a singular 
subarc at the end of the battle, it cannot be an optimal 
trajectory. 

From studying this problem we feel that following result 
is true. . If so it would reduce this particular problem to 
the two cases for the optimum, i.e., u = 0 or u h 1. We 
were unable to prove this result or to find a counter example 
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THEOREM 5 



A maximizing arc never has a singular subarc nor a corner 
that is the optimal control is always either u=lorusO, 
depending on the given conditions. 

(Not proved.) 



\ 
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VI. DISCUSSION 



Considering the theorems given above, we see that if an 
extremal does not contain a singular subarc, the optimal 
control u* is (a) constant for all te(0,T) and is (b) either 
0 or 1. Which of these policies is better may be determined 
by direct computation and comparison of values of the cri- 
terion functional using the attrition-rate coefficients and 
the initial conditions. Thus, we see that in this case the 
optimal time-sequential fire-support policy for Blue is to 
concentrate all his artillery resources on one of the Red 
units for the entire duration of battle. It should be 
pointed out that even if an extremal contains a singular 
subarc. Blue's artillery fire will (optimally) still be 
concentrated on one of the Red units at the end of battle. 
This is because an optimal fire-support policy has been 
characterized as not terminating with a singular subarc. 

We have not been able to show which policy is optimal 
when the trajectory contains a singular subarc for some 
period of time before the end of the problem. Thus, the 
complete explicit determination of the optimal fire support 
policy was not accomplished. However, the basis has been 
given for doing this in the future. 

It should be pointed out that the insights gained into 
optimal fire-support policies obtained from the model (1) 
are no more valid than the model itself. Thus the reader 
should be cautioned that the optimal fire-support policy 
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determined here may only apply to the problem (1). As 
discussed by Taylor (see [14]), a combat optimization prob- 
lem consists of the following three parts: (1) decision 

criteria, (2) model of the planning horizon (e.g. conflict 
termination condition), and (3) model of combat dynamics. 
For different objectives or different combat dynamics, the 
optimal fire -support policy may well be different from that 
given here. 
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VII. CONCLUSIONS 



For the problem considered in this thesis, the optimal 
fire-support policy is piecewise constant with at most two 
switches. If one never divides fire between the two enemy 
units (singular control) , then the optimal policy is con- 
stant and either 0 or 1 (i.e., concentrate all fire on one 
of the enemy units) . Although complete details could not 
be worked out within the time constraints for this thesis, 
it is conjectured that this is the optimal policy (i.e., 
always concentrate all supporting fire on one of the enemy 
units). The reader should be cautioned that these conclu- 
sions may only apply to the specific fire support problem 
studied here. 
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